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Abstract. 

Recent results from BES and CLEOc experiments on hadron spectroscopy and charmonium 
decays using i^' and \(/" data samples collected in e+e^ annihilation are reviewed, including 
the observation of X(1835) in 7/1/ ^ /tt+tt^T]', study of the scalar particles in J/xj/ radiative and 
hadronic decays, as well as in XcO hadronic decays, and the study of the "pTT puzzle" in J/y/, xj/', 
and decays. 
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1. INTRODUCTION 

BESII [1] running at BEPC and CLEOc [2] running at CESR are the two detectors 
operating in the T-charm energy region, and have collected large data samples of char- 
monium decays including 58 M 7/i// events, 14 M x^f' events, and 33 pb^^ data around 
1//" peak at BESII, and 4 M i//' events, and 281 pb^^ i//" events at CLEOc. To study 
the continuum background in the charmonium decays, special data samples at center of 
mass energy lower than the i//' mass were taken both at BESII {y/s = 3.65 GeV) and at 
CLEOc {^/s = 3.671 GeV), the luminosity are 6.4 pb^^ and 21 pb^^ respectively. These 
data samples are used for the study of the hadron spectroscopy, the D decay properties 
and the CKM matrix, as well as the charmonium decay dynamics. 

In this paper, we focus on the search for the new resonances in 7/i// decays, the 
properties of the scalars from 7/ 1// radiative and hadronic decays, and a new approach of 
studying the scalars using Xcj decays, and the extensive study of the "p;r puzzle" related 
physics in 7/ 1//, i//' and x^r" decays. 

It should be noted that the CLEOc detector is much better than the BESII detector, 
especially in the photon detection, this makes its 4 M i//' events data sample produces 
results with similar precision as from 14 M i//' events from BESII. 



2. OBSERVATION OF X( 1835) 

The decay channel 7/1// yn^n^rj', with rj' k^k^t], or r\' yp, is analyzed 
using a sample of 58 x 10^ 7/1// events collected with the BESII detector [3], to search 
for the other decay modes of the possible existing pp bound state as observed in 
J /\lf ^ ypp process at BESII [4]. Figure 1 shows the tz^tz^T]' invariant mass spectrum 
for the combined 7/i// yn^ r\' {r\' n^n^T]) and 7/i// yn^ r\' {r\' yp) 
samples. A clear peak is observed at around 1 .835 GeV/c^. The spectrum is fitted with a 
Breit-Wigner (BW) function convolved with a Gaussian mass resolution function (with 
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(7=13 MeV/c^) to represent the X(1835) signal plus a smooth polynomial background 
function. The signal yield from the fit is 264 ± 54 events and the statistical significance 
is 7.7 o. 

120 I 1 




M(7tVr|') (GeV/c^) 

FIGURE 1. The n^Tt^rj' invariantmass distribution for selected 7/ y/^ Yn^n~ri'{ri' ;r+;r"77, 77 
YY) and J/y/ ^ y^^^ V'iv' ~^ YP) events. Dots with error bars are data, sohd line is the fit, and the 
dashed curve indicates the background function. 

The mass of X(1835) is measured to be M = \S33 J ±6.l{stat) ±2.7 (syst) MeV/c^, 
the width is F = 67.7 ±20. 3(stat) ±1.1 (syst) MeV/c^, and the product branching 
fraction is 5(7/1//^ 7X) ^ Tt+Tt-T]') = {2.2±0.4{stat) ±0.4{syst)) x IQ-^^. The 
mass and width of the X(1835) are not compatible with any known meson resonance 
listed by PDG [5]. In Ref. [5], the candidate closest in mass to the Z(1835) is the 
(unconfirmed) 2"+ T]2 ( 1 870) with M = 1 842 ± 8 MeV/c^. However, the width of 225 ± 
14 MeV/c^, is considerably larger than that of the X(1835) (see also [6], where the 2 ^ 
component in the 77 TtTt mode ofJ/\j/ radiative decay has a mass 1840±15MeV/c2anda 
width 170 ± 40 MeV/c^). Another candidate with mass close to the X ( 1 835) in PDG [5] 
is the Z(1860) observed in the pp mass threshold in radiative J/^f ^ ypp decays [4], 
where a mass of 1859^27 MeV/c^, and width smaller than 30 MeV/c^ at 90% C.L. were 
reported. It has been pointed out that the S-wave BW function used for the fit in Ref. [4] 
should be modified to include the effect of final- state-interactions (FSI) on the shape of 
the pp mass spectrum [7, 8]. By redoing the 5-wave BW fit to the pp invariant mass 
spectrum of Ref. [4] including the Isospin zero, S-wave FSI factor of Ref. [8], BES 
gets a mass M=1831±7 MeV/c^ (in good agreement with the X(1835)) and a width 
r < 153 MeV/c^ at the 90% C.L. (not in contradiction with the X(1835)). There are also 
theoretical arguments that the two states, Z( 1835) and X( 1860), are indeed one particle, 
and is a pp bound state [9, 10], however, other possible interpretations of the Z(1835) 
that have no relation to the p'p mass threshold enhancement are not excluded. Further 
information about the X(1835) and X(1860), especially better precisions on the mass 
and width measurements and the determination of the spin-parity are strongly desired 
before concluding the nature of the states. 
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3. SCALAR PARTICLES IN J/\\f AND Xco DECAYS 

The study of the scalars are very important in two aspects: in experiment, there are still 
controversies about the resonance parameters of these states; and in theory, it is still hard 
to incorporate all the experimental results in a self-consistent picture. The reason for the 
former is somewhat due to the techniques used in extracting the physics information 
from the experimental data, namely, the partial wave analysis (PWA) was extensively 
used in the analyses, but sometimes it is rather arbitrary what resonance need to be 
included in the complicated fit with a few ten and even more than one hundred free 
parameters. The reason for the latter, in part is due to the fact that the experimental 
results may not all be reliable, and the complexity in the low energy QCD domain 
that the higher order terms neglected may not be small, and the mixing of the states 
in principle is hard to be considered completely. 

3.1. Radiative and hadronic J/x^f decays 

Using the world largest //i// data sample in e^e~ annihilation experiment, BES 
studied the scalars decay into pair of pseudoscalars {n^n^ , n'^n^, K^K and K^K^) 
inJ/xi/ radiative decays as well as recoiling against a or an ft) [1 1, 12, 13, 14]. The full 
mass spectra and the scalar part in them are shown in Fig. 2. 

From the analyses, BES sees significant contributions of a particle in (on^n'^ and 
(oK^K , and also hint in ^k^k^. Two independent partial wave analyses are performed 
on (OTf^Tf^ data and four different parameterizations of the o amplitude are tried, all give 
consistent results for the o pole, which is at (541 ± 39) - /(252 ± 42) MeV/c^. There 
is also events accumulation in the low 7l^n~ mass in Y^~^^~ mode, most probably 
due to the contribution of the a, however, there is no attempt to analyze the structure 
at BES, one possible reason is the presentation of the huge background from p^Tt^. 
Nevertheless, the coupling of the a with a photon in 7/ Xj/ decays is an important piece 
of information for the understanding of the nature of the particle, a better detector with 
more statistics may want to measure it. A measurement of the resonance parameters 
of G using the BESII i/ — > J/\i/7t'^7t~ data gives similar pole position as measured in 
J/Y^ (on+n~ [15]. 

Strong /o(980) is seen in ^K^n^ and ^K^K modes, from which the resonance pa- 
rameters are measured to be M = 965 ±S{stat) ±6{syst) MeV/c^, gi — 165 ± I0{stat) ± 
I5{syst) MeV/c^ and gi/gi = 4.21 ±0. 25 {st at) ±0.21(5');5'?), where Mis the mass, and 
gi and g2 are the couplings to nn and KK respectively if the /o(980) is parameterized 
using the the Flatte's formula. The production of /o(980) is very weak recoiling against 
an ft) or a photon, which indicates ss is the dominant component in it. 

The <I>'K'^'K~ data also show a strong scalar contribution at around 1.4 GeV/c^, it 
is due to the dominant /o(1370) interfering with a smaller /o(1500) component. The 
mass and width of /o(1370) are determined to be: M = 1350 ±50 MeV/c^ and T = 
265 ±40 MeV/c^. In Y^^^~, a similar structure is observed in the same mass region, 
the fit yields a resonance at mass 1466 ±6{stat) ± I6{syst) MeV/c^ with width of 
108t\ilst at) ±21 (sy St) MeV/c^, possibly the /o(1500), and the contribution from the 
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FIGURE 2. The invariant mass distributions of the pseudoscalar meson pairs recoiling against co, (p, or 
yinJ/y/ decays measured at BESII. The dots with error bars are data, the solid histograms are the scalar 
contribution from PWA, and the dashed lines in (a) through (c) are contributions of a from the fits, while 
the dashed line in (d) is the /o(980). Notice that not the full mass spectra are analyzed in (e), (f), and (g). 
Results in (e) are preliminary, otherwise are published. 



/o(1370) can not be excluded. The production of /o(1370) and /o(1500) in jKK is 
insignificant. 

The K^K invariant mass distributions from yKK and coK^K , the n^n^ in- 
variant mass distributions from yn^TT^, and ^tt^tt^ show clear scalar contribution 
around 1.75 GeV/c^. Two states are resolved from the bump, one is /o(1710) with 
M ~ 1740 MeV/c^ and T ~ 150 MeV/c^ which decays to KK mostly, and one possi- 
ble new state /o(1790) with M ~ 1790 MeV/c^ and T ~ 270 MeV/c^ which couples 
to TTTT Stronger than to KK. However, the existence of the second scalar particle needs 
confirmation: the signal observed in 0/o(179O) is rather in the edge of the phase space. 
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and the reconstruction efficiency of the decreases dramatically as the momentum of 
the decreases thus the momentum of the kaon from ^ decays is very low and can not 
be detected. Furthermore, there are wide higher mass scalar states above 2 GeV/c^ as 
observed in 7;r+;r~ (Fig. 2e) and yKK [5], whose tails may interfere with the /o(1710) 
and produce structure near the edge of the phase space. 

The use of these measurements for understanding the nature of the scalar particles can 
be found in Refs. [16, 17, 18, 19], where the J/ x^f decay dynamics and the fractions of 
the possible qq and glueball components in the states are examined. 

3.2. Pair production of scalars in XcO 7i^7i~K^K~ 

Partial wave analysis of XcQ ^ k^k K+K is performed [20] using XcQ produced in 
x^f' decays at BESII. In 14 M produced x^f' events, 1371 i//' yxm, XcO n^n~K^K~ 
candidates are selected with around 3% background contamination. 

Fig. 3(a) shows the scatter plot of K^K" versus ti'^'K~ invariant mass which provides 
further information on the intermediate resonant decay modes for {n'^n~){K'^K~) 
decay, while Fig. 3(b) shows the scatter plot of K^K~ versus K~K^ invariant masses 
for the resonances with strange quark. 

Besides {'K'k){KK) and {K'k){K'k) modes, {KTnt)K mode which leads to a measure- 
ment of ^1 (1270)^ and (1400)^ decay processes is also included in the fit. The PWA 
results are summarized in Table 1. From these results, we notice that scalar resonances 
have larger decay fractions compared to those of tensors, and such decays provide a 
relatively clean laboratory to study the properties of scalars, such as /o(980), /o(1370), 
/o(1710), and so forth. The upper limits of the pair production of the scalar mesons 
which are less significant are determined at the 90% C.L. to be 
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FIGURE 3. The scatter plots of (a) K^K~ versus TT+TT" and (b) K'^7t~ versus K~n^ invariant mass 
for selected xj/' yXcO, XcO ^ 7C+n~K+K~ events. 
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^/o(1500)/o(1370)]=^[/o(1500) ^ K+K-\m[fomQ) ^ K+R-] < 1.4 x IQ-^, 
■^[Zco ^/o(l 500)/o( 1500)]^ [/o( 1500) ^ 7r+7r"]^[/o(1500) ^ K+R-] < 0.55 x 10"^, 
^\XcO /o(1500)/o(1710)]^[/o(1500) ^ ;r+;r-]^[/o(1710) ^ K+K \ < 0.73 x lO"'*. 



TABLE 1. Summary of the XcO n^n^K^K'^ results, where X represents the 
intermediate decay modes, A^-^" is the number of fitted events, and s.s. indicates signal 
significance. 



Decay mode 




Branching ratio (10 "*) 


S.S. 


(X) 




^[Xco -^X-f n+n-K+R-] 




/o(980)/o(980) 


27.9 ±8.7 


3.46 ±1.0811:^^ 


5.3a 


/o(980)/o(2200) 


77.1 ±13.0 


8.42±1.42t[f^ 


7.1a 


/o(1370)/o(1710) 


60.6 ±15.7 


7.12±1.85;^-2^ 


6.5a 


(892)0^* (892)" 


64.5 ±13.5 


8.09 ±1. 69;^ 


7.1a 


A'(J(1430)£,](1430) 


82.9 ±18.8 


10.44 ±2.37lf^^ 


7.2a 


A'i;(1430)^^(1430) +C.C. 


62.0±12.1 


8.49±1.66t||^ 


8.7a 


Ki{\21Q)+K- +C.C., 




9.32±1.83l};^i 




Ki{\210) ^ KpillO) 


68.3 ±13.4 


8.6a 


Kx{UOO)+K- +C.C., 








Ki{U00)^K*{%92)n 


19.7 ±8.9 


< 11.9 (90%C.L.) 


2.7a 



The above results supply important information on the understanding of the natures 
of the scalar states [21]. 



3.3. Comments on the PWA 

PWA is extensively used in extracting physics information from the experimental data, 

all the information listed above in this section is from PWA. While we know the PWA 
uses the information in the full phase space for physics study so that it is more powerful 
than working in one dimension (invariant mass, for example) or a bit higher dimension, 
it suffers from too many free parameters and other technical problems. 

First of all, the experimental data are all affected by the finite resolution in momentum, 
energy, or direction measurement, this was not considered in current analyses; secondly, 
the parametrization of the resonance, especially the wide resonance, still have room to 
improve. Finally, the effect of the imperfect simulation of the detectors is hard to be 
considered in a fit with a few ten or even more free parameters. These effects may not 
be significant when the statistics is low, however, as the statistics increases, all these 
effects will possibly produce significant fake signals. How to handle these, if not now, 
at least in the near future, should be studied since high luminosity experiments will be 
soon available. 

Studying the papers dealing with the PWA, it is found that two important information 
are missing in most of the analyses, which are the goodness-of-fit and the correlation 
coefficients between the fit out parameters. 

In many of the analyses, the comparison to the alternative fits are given to show the 
fit is the best among all the fits tested, however this does not guarantee the fit used in the 
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analysis is reliable, just like a fit to a Gaussian using a 2nd order polynomial is better than 
using a straight line, but does not mean the fit is acceptable. To give the goodness-of-fit is 
not easy in case of PWA since the fitting function is always very complicated. A possible 
way is to supply a simple j}' test in one or multi-dimension as has been done in Ref. [20], 
although not perfect, it shows the reader a feeling how the fit describes the data. Another 
possible way is try to use toy Monte Carlo to get the expected distribution of the or the 
likelihood after generating many MC samples using the fit out parameters as input, and 
compare with the one got from the fit to the data — this may be a bit time consuming, 
however, it is more reliable since the detector effects are considered. 

The parameters from the fit are generally correlated, and sometimes, some variables 
are highly correlated, in this case, only reporting the diagonal error is not enough. This 
is extremely important when there are two nearby resonances with strong interference, 
they are anti-correlated and one component will increase (decrease) as the other decrease 
(increase) to make the total contribution unchanged. Neglecting this in the theoretical 
analysis will also be dangerous. Another effect of the correlation is the significance of 
the signal may be affected strongly, that is, a declared significance may not be as high 
as that determined when there is no correlation, and vice versa. 



4. "p;r PUZZLE" IN VECTOR CHARMONIA DECAY 

From perturbative QCD (pQCD), it is expected that both J and ^' decaying into light 
hadrons are dominated by the annihilation of cc into three gluons or one virtual photon, 
with a width proportional to the square of the wave function at the origin [22]. This 
yields the pQCD "12% rule", 

Qh = = f"^"^^"^' ^ 12%. (1) 

A large violation of this rule was first observed in decays to pn and K*^K +c.c. 
by Mark II [23], known as the pit puzzle, where only upper limits on the branching 
fractions were reported in i/a' decays. Since then, many two-body decay modes of the I// 
have been measured by the BBS collaboration and recently by the CLEO collaboration; 
some decays obey the rule while others violate it [24, 25]. 

The extension of the above rule to i/' is straightforward, in the same scheme, one 
would expect 

Q!h = §^ = y-'"'- ^ 1.9 X 10-4. (2) 

As the i/a" data samples are available both at BESII and CLEOc, the search for the 
decays of y/" into light hadrons was performed. Since y/" is above the charm threshold, it 
is expected the dominant decays of it is to charmed meson pairs, however, large fraction 
of charmless decays of \]/' is expected if x/' is a mixture of S- and Z)-wave charmonium 
states and the mixing is responsible for the "12% rule" violation in J/ \\r and i/a' decays. 
The above two rules may be tested by the large data samples ofJ/xj/, V ' ^^^^ V ^t both 
BESII and CLEOc. 
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4.1. \\f ^ pjz: current status 

The p% mode of the vector charmonia decays is essential for this study, since this is 
the first puzzling channel found in 7/ 1//^ and i// decays. The new measurements, together 
with the old information, show us a new picture of the charmonium decay dynamics. 



4.1.1. J/^f^Tf^n 7C^ 

BESII measured the 7/1// — » k^k^k'^ branching fraction using its 7/i/a and i/a' 
J / \i/n^n^ data samples [26], and BARBAR measured the same branching fraction us- 
ing 7/ 1// events produced by initial state radiative (ISR) events at = 10.58 GeV [27]. 
Together with the measurement from Mark-II [23], we get a weighted average of the 
^{J/\l/^n+n-7t^) = (2.00±0.09)%. 

To extract the J/y/ ^ p7t branching fraction, PWA is needed to consider the possible 
contribution from the excited p states, the only available information on the fraction of 
pit in J/ Y ^ Tl^ll^ll^ was got at Mark-IE. Using the information given in Ref. [28], 
we estimate ^^'^^^J^^-^^o-) = 1-1'7(1 ± 10%), with the error from an educated guess 
based on the information in the paper since we have no access to the covariant matrix 
from the fit showed in the paper. From this number and the SS{J j \j/ — > n^n^n^) got 
above, we estimate ^{J/y ^ P^r) = (2.34 ±0.26)%. This is substantially larger than 
the world average Usted by PDG [5], which is ( 1 .27 ± 0.09) %. 



4.1.2. Iff' ^ 71+71-71^ 

X]/ ^pTt was studied both at BESII [29] and CLEOc [25] . After selecting two charged 
pions and two photons, clear tt^ signals are observed in the two photon invariant mass 
spectra, the numbers of signals are found to be 229 and 196 from BESII (shown in left 
plot of Fig. 4 as an example) and CLEOc respectively, and the branching fraction of 
Xj/'^n+n'TC^ is measured to be (18.1 ± 1.8± 1.9) x 10-^ and (18.8+} J ± 1.9) x 10-^ 
at BESII and CLEOc respectively. The two experiments give results in good agreement 
with each other. The Dalitz plot of y' K^7U~7l^ events (right plot of Fig. 4) shows 
very different signature from that of J /y Tf^Tt'Tt^ decays, there are lots of events in 
the center part of the Dalitz plot in the former case, while in the latter, almost all the 
events are in the p mass region and there is almost nothing in the center of the DaUtz 
plot. 

To extract the branching fraction of X]/ ^ pn, however, BESII uses a PWA including 
the high mass p states and the interference between them, while CLEOc counts the 
number of events by applying a p mass cut. The branching fraction from BESII is 
(5.1 ±0.7 ±1.1) X 10"^ while that from CLEOc is (2.4+[];5±0.2) x 10"^ the difference 
is large. Although big difference exist between BESII and CLEOc results exist, these do 
mean the i// — > p;r signal exists, rather than the signal is completely missing as has been 
guessed before. If we take a weighted average neglecting the difference between the two 
measurements, we get —>■ pn) = (3.1 ±0.7) x 10~^. 
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FIGURE 4. Two photon invariant mass distribution (left) and the DaUtz plot (right) after final selection 
for BESII data. The histograms are data, and the curves show the best fits. 



Comparing ^ pn) with 3§{J/\j/^ pn), one gets 



Qpn 



The suppression compared to the 12% rule is obvious 



(0.13 ±0.03)%. 



4.1.3. Xj/" ^7C+7C-7C^ 

It has been pointed out that the continuum amplitude plays an important role in 
measuring decays into light hadrons [30]. In fact, there are two issues need to be 
clarified in decays, that is whether decays into light hadrons really exist, and if 
it exists, how large is it. By comparing the cross sections of e^e~ n~^7t~7t^ at the i/a" 
resonance peak (^/s = 3.773 GeV) and at a continuum energy point {^/s — 3.65 GeV 
at BESII and 3.671 GeV at CLEOc) below the peak, both BESH and CLEOc found 
that <7(e+e~ — > 7t'^7t~7t^) at continuum is larger than that at y/" resonance peak. The 
average of the two experiments [31, 32] are 

(7(e+e~ ;r+;r~;r°)on = 7.5±1.2pb, 

a{e+e~ 7t+7t~n^)Qff = 13.7 ± 2.6 pb. 

The difference, after considering the form factor variation between 3.65 and 3.773 GeV, 
is still significant, and it indicates there is an amplitude from i//" decays which interferes 
destructively with the continuum amplitude, which makes the cross section at the y/" 
peak smaller than the pure contribution of continuum process. 

For the pn mode, BESII can only give upper limit of its cross section due to the lim- 
ited statistics of the data sample, the upper limit at 90% C. L. is found to be 6.0 pb [31] 
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at the y" peak, which is in consistent with the measurement from CLEOc using a much 
larger data sample: a(e+e^ p^)on = 4.4±0.6 pb [32]; while the cross section at the 
continuum is 8.0_i'4 ± 0.9 pb measured by CLEOc. 

To extract the information on the Xjr" ^ pn branching fraction, BESII developed a 
method based on the measured cross sections at y/" resonance peak and at the con- 
tinuum [31]. By neglecting the electromagnetic decay amplitude of i//", there are two 
amplitudes contribute to the cross section at the i//" peak, the strong decay amplitude 
of y/" and the continuum amplitude. Taking the continuum amplitude as a real number, 
the Xj/' strong decay amplitude is described as one real number for the magnitude, and 
one phase between i//" strong and electromagnetic decays to describe the relative phase 
between the two amplitudes. Since only two measurements are available (at yr" peak 
and at continuum), one can only extract y" decay branching fraction as a function of 
the relative phase. BESII measurement on the upper limit of the e+e ^ pTZ cross sec- 
tion at yf" peak, together with the CLEOc measurement of the continuum cross section 
restrict the physics region of the branching fraction and the relative phase as shown in 
Fig. 5(left). From the Figure, we see that the branching fraction of yr" ^ pTZ is restricted 
within 6 x 10"^ and 2.4 x 10"^ and the phase is between -150° and -20°, at 90% C.L. 




-ISO -160 -140 -120 -100 -80 -60 -40 -20 



Phase between Vf" Strong and EM decays 

FIGURE 5. Physics region on pn) and the relative phase between y/" strong and electromag- 

netic decays from BESII (left); and the illustration of the two solutions in y/" decays (right): OA represents 
the continuum amplitudes, OB or OC represents the peak amplitudes, and AB and AC are the two solutions 
for the resonance decay amplitudes, a is the relative phase between continuum amplitude and the y/" 
strong decay amplitude, including the relative phase and the relative phase between continuum and y/" 
Breit-Wigner amplitudes. 

The observation of the e+e^ pK signal at yf" peak and the measurement of the 
cross section [32] at CLEOc further make the physical region in the branching ratio 
and relative phase plane smaller: the CLEOc measurement gives a similar out bound of 
the physical region as BES gives, while also indicates the central part of the physical 
region in Fig. 5(left) is not allowed by physics. By using a toy Monte Carlo to simulate 
the CLEOc selection criteria and the interference between resonance and continuum 
amplitudes, we found that the yf" ^ pn branching fraction could be either (2. 1 ± 0.3) x 
10^^ or (2.4^2 o) ^ 10^^ from the CLEOc measurements, if the relative phase between 
yr" strong and electromagnetic decay amplitudes is —90° as observed in J /yr and yr' 
decays [33]. 
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The reason why there are two solutions for x^f" decays can be understood as illustrated 
in Fig. 5b. If we take the continuum amplitude as a real number, it can be shown as 
a vector along the real axis in the complex plane, the total cross section at y" P^ak 
only gives the magnitude of the total amplitude and it is shown as a circle in the plot, 
the amplitude of the i//" decays, represented by a vector connecting the end of the 
continuum amplitude and the end of the total amplitude may have two cross points 
with the circle, representing the two solutions of the l//' decay amplitudes, and thus 
the branching fractions. Only in some very special cases, there is only one solution. The 
angle between the continuum amplitude and the decay amplitudes is shown in the 
plot as a, it is the sum of the relative phase ^ between the strong and electromagnetic 
decays of Xjf", and the phase from Breit-Wigner function for the resonance. It can be 
seen that, the two solution only happens when the cross section at continuum is larger 
than or equal to that at i//" peak; otherwise, there is only one solution, as in x^f' ^ pn 
case. However, as in physics, there is only one Xj/" decay branching fraction, there must 
be a way to extract the ij/" decay branching fraction experimentally, this could be made 
possible if one more data sample is taken at a different energy point, for example, in 
the shoulder of the i//" resonance; furthermore, if one even wants to get the value of 
the relative phase, one more data point is necessary. In total, to get concrete physics 
information on the i/' decay branching fraction, at least data samples at four different 
energy points are needed, better have one point far from the resonance, so that it gives 
absolute measurement of the continuum amplitude. 

Based on current data samples, one gets 

to be compared with the pQCD prediction of 0.019% if i/' is pure D-wave charmonium. 



4.2. Other studies and comments 

There are many more new measurements on \]/ decays for the extensive study of 
the "12% rule" [24, 34, 25], among which the Vector-Pseudoscalar (VP) modes are 
measured at first priority. The ratios of the branching fractions are all suppressed for 
these VP modes compared with the 12% rule. The multi-hadron modes and the Baryon- 
antibaryon modes are either suppressed, or enhanced, or normal, which are very hard 
to categorize simply. The various models, developed for interpreting specific mode may 
hard to find solution for all these newly observed modes. 

The y/" decays into light hadrons were searched for in various y/" decay modes, 
including VP and multi-hadron modes [32, 35]. However, only the comparison between 
the cross sections at continuum and those at resonance peak are given, instead of 
giving the yr" decay branching fractions. In current circumstances, it is still not clear 
whether the y/" decays into light hadrons with large branching fractions, since, as has 
been shown in the pit case, there could be two solutions for the branching fraction, and 
the two values could be very different. 
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As the luminosity at the y" peak is large enough, current study is limited by the low 
statistics at the continuum: at CLEOc, the luminosity at continuum is more than an order 
of magnitude smaller than that at peak, this prevents from a high precision comparison 
between the cross sections at the two energy points. One conclusion we can draw from 
the existing data is that the measurements do not contradict with the assumption that 
the relative phase between y/" strong and electromagnetic decay amplitudes is around 
—90°, and the X]/' decays into light hadrons could be large. 

The study of the pn puzzle between i/a' and J/yf decays and the charmless decays of 
should not be isolated as they were since //v^, V^' and i/a" are all charmonium states 
with very similar quantum numbers, and it is expected and are the mixtures of 2S 
and ID states [36]. In developing models to solve one of the problems, the others should 
also be considered. There have been a few models developed following this line or can 
be easily extended to all these three states, like the S- and D-wave charmonia mixing 
model [36], the DD re- annihilation in [37], the four-quark component in Xj/" [38], 
and survival cc in \]/' [39], and so on. Experimentally testable predictions are welcome 
for validating the models. 

One further observation is that many of the attempts to interpret the pn puzzle are 
based on the potential models for the charmonium which were developed more than 20 
years ago, as the B-factories discovered many new charmonium states [40] which are 
hard to be explained in the potential models, it may indicate even the properties ofJ/\j/, 
and y/" are not as expected from the potential models. The further understanding of 
the other high mass charmonium states may shed light on the understanding of the low 
lying ones. 

5. SUMMARY 

There are many new results on hadron spectroscopy from the charmonium decays from 
BBS and CLEOc experiments. While many analyses supply more information on the 
known states like the light scalar particles to understand the nature of them, there are 
also new observations which may indicate there are still something unexpected in the 
low energy regime such as the possible existing baryonium. The decay properties of the 
vector charmonium states, although have been studied for more than three decades, is 
still far from being understood, one extreme example is the "p;r puzzle" inJ/yr and yf' 
decays. Further studies of all these are expected from the BESni at BEPCn which will 
start its data taking in 2007. 
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